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Resumen 

El cerro Pico de Piedra constituye una estructura volcánica colapsada del Grupo volcánico Irazú, 
asociable con el vulcanismo del miembro Paraíso y conforma la cuenca alta del río Virilla. La 
interpretación por medio de sensores remotos, recopilación bibliográfica, levantamiento geomorfo- 
lógico y visitas al campo permitieron identificar procesos erosivos intensos en la cuenca, inestabi- 
lidad de laderas y susceptibilidad del corredor vial La Holanda que comunica fincas productoras de 
leche y tubérculos con San Isidro de Coronado. La ocurrencia de deslizamientos y susceptibilidad 
se evaluó con el método Mora-Vahrson modificado y Probabilistic Infinite Slope Analysis (PISA), 
modelando diversos grados de saturación en los suelos volcánicos y diversos escenarios, entre 
ellos los de Patillos de 1952 y Cinchona en 2009. La zonificación geológica integrada en la cuenca 
sugiere conservar la cobertura forestal actual, controlar la erosión y respetar los retiros respecto a 
los cauces fluviales. Las áreas agropecuarias actuales son de uso conforme y las más estables en 
relación a su condición geomorfológica. 


Palabras clave: Deslizamientos, suelos volcánicos, arrastre de bloques, amenaza sísmica, 
geomorfología. 


Abstract 
The Cerro de Pico de Piedra (translated as Mount Pico de Piedra) is a collapsed volcanic structure 
of the Irazú Volcanic Group associated with Paraíso Member”s volcanism and constitutes part of 
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the upper basin of the Virilla River. Data interpretation gathered by remote sensing, bibliographical 
compilation, geomorphological survey and field visits allowed the identification of intensive ero- 
sion processes, slope instability and hazard (damage) susceptibility ofthe La Holanda transport co- 
rridor (road) that connects milk- and tuber-producing farms with the town of San Isidro de Corona- 
do. Landslide occurrence and transport corridor susceptibility were assessed according to modified 
Mora-Vahrson and probabilistic infinite slope analysis (PISA) methodologies by modeling various 
degrees of soil water saturation in volcanic soils under two conditional seismic scenarios referenced 
to the Patillos (1952) and Cinchona (2009) Costa Rican earthquakes. Integrated geological zoning 
of the upper basin suggests the need to preserve the current forest coverage, minimize erosion rates 
and maintain appropriate buffering allowance distances of anthropogenic activities from fluvial 
bodies. Current agricultural areas are located within the legally established regulated territorial 
zoning attributed as being the most stable in relation to their geomorphological conditions. 


Keywords: Landslides; Volcanic soils; Boulder drag; Seismic hazard; Geomorphology. 


Introducción 

El área de estudio se localiza en el cantón Vázquez de Coronado, 
provincia de San José. Abarca el sector este del distrito de San Rafael, 
siendo un sector donde se asientan fincas de producción agropecuaria, in- 
fraestructura hidráulica de captación y distribución, además se localizan 
algunas viviendas (Figura 1). En este sector inicia la cuenca del río Virilla 
y de la quebrada Sangrino, de las cuales se capta agua para algunos pobla- 
dos del cantón. 

Predominan las laderas abruptas, sectores boscosos y pastos dedi- 
cados a producción pecuaria y en menor medida a parcelas de cultivo. El 
punto más elevado es el cerro Pico de Piedra con 2633 m.s.n.m. y consti- 
tuye una estructura volcánica colapsada asociada al volcanismo del Grupo 
Irazú (Alvarado et al., 2006, p. 263). 

La precipitación anual alcanza los 2500 mm, siendo septiembre 
y octubre los meses más lluviosos (IMN, 2019), además se presentan 
suelos descubiertos y evidencias superficiales de inestabilidad en suelos 
volcánicos. 

Los andisoles relacionados al Volcán Irazú (Mata, Rosales, Vás- 
quez & Sandoval, 2013) presentan texturas franco arenosas a arcillosas, 
con buena retención de humedad, baja densidad aparente y baja resis- 
tencia al corte o cizalle (Henríquez, Cabalceta, Bertsch & Alvarado, 
2014). La fábrica de estos suelos y el intemperismo de materiales volcá- 
nicos recientes usualmente genera alófanas y halloysitas, las cuales son 
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Figura 1: Uso de la tierra dentro del área de estudio en las cercanías del 
cerro Pico de Piedra en Vázquez de Coronado, los límites del área y su 
localización dentro de Costa Rica. 
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Fuente: Insumos cartográficos IGN (2007) y Sistema Nacional de Información Territorial (2017) 


susceptibles ante la erosión hídrica acelerada (Alvarado, Mata  Chinchi- 
lla, 2014, p. 82 y 86). 

Por otra parte, el acceso al sitio se logra por medio de la red vial can- 
tonal (en adelante corredor vial La Holanda), la cual cuenta con un puente 
de propiedad privada, dos puentes viales cantonales sobre el cauce del río 
Virilla (P1 y P2), tendidos eléctricos y de telecomunicaciones (Figura 2) 
susceptibles ante colapsos por procesos de remoción en masa y caída de 
rocas, principalmente. 
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Figura 2: Mapa inventario de infraestructura en el área de estudio 
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Adicionalmente, el inventario de infraestructura indica una calzada 
vial con un ancho de 4 m (líneas negras y grises Figura 2), cinco tanques 
de captación y sus tuberías pertenecientes al Instituto Costarricense de 
Acueductos y Alcantarillados (AyA), y más de cincuenta inmuebles entre 
viviendas, almacenes, graneros y lecherías. 


Contexto geológico y sísmico 

En primera instancia, Bergoeing (1979, p. 139 y 140) identifica una 
estructura volcánica que denominó volcán Las Nubes, que consiste de un 
antiguo cráter erosionado en el extremo oeste (Figura 1). 

Vargas (1994, p. 37 y 38) estudia el área con fines hidrogeológicos 
determinando que las rocas aflorantes en el área de estudio son basaltos 
y andesitas con augita asociados con el Grupo Aguacate. Sin embargo, 
investigaciones como Alvarado (1993) y Alvarado et al. (2006, p. 264) 
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describen la Unidad Pico de Piedra con una edad de 569 ka (Chibania- 
no) y que constituye las rocas más antiguas sobre las que evolucionó el 
estratovolcán actual Irazú. Se trata de andesitas y andesitas basálticas 
con alteración hidrotermal severa que afloran en valles fluviales, algunos 
puntos de las colinas y en el cerro Pico de Piedra, donde la cobertura de 
suelos es escasa. 

La Unidad Pico de Piedra es considerada como parte de la Formación 
Reventado (etapa temprana del volcanismo del Irazú) y a su vez dentro 
de todos los depósitos, piroclastos y lavas del Grupo Irazú (Krushensky, 
1972), que cubren el sector norte de la provincia de Cartago, así como el 
margen sureste y este del cantón de Vázquez de Coronado (Figura 3). 

Sobreyaciendo estas rocas y depósitos se identifican cenizas recien- 
tes y tobas sin diferenciar que se extienden desde el foco eruptivo del Irazú 
alcanzando localidades como Platanares de Moravia, San Isidro, Cascajal, 
y Las Nubes de Coronado (Vargas, 1994, p. 43; Denyer & Arias, 199la, 
p. 45 y 46). La alteración de estos materiales ha generado la cobertura de 
suelo actual. 

Otras hipótesis como la planteada por Hidalgo, Linkimer & Alva- 
rado (2004, p. 106), determinan que la abertura al oeste de la estructura 
volcánica Pico de Piedra y la estructura Cabeza de Vaca se asocian con un 
evento o serie de eventos de colapso gravitacional, cuyos depósitos alcan- 
zaron el Valle Central como avalanchas (debris avalanches/debris flows) 
de rocas y partículas finas, cargadas en agua y constituyeron el “aglutinan- 
te” de los materiales, expuestos en la mayoría de cauces fluviales del Valle 
Central (Figura 3). 

Se identifican otras rocas en el margen oeste de la Figura 3 que 
forman parte de la Paleo-Cordillera, la cual incluye los antiguos comple- 
jos volcánicos del Grupo Volcánico Zurquí (Obando, 2004, p. 84 y 85; 
Kussmaul, 2000, p. 83 y 84), anteriores a los volcanes actuales de la Cor- 
dillera Volcánica Central. En el margen este por su parte, afloran rocas 
contemporáneas a las del área de estudio asociadas al volcán Turrialba, 
entre ellas lavas y piroclastos (Alvarado et al., 2006, p. 263; Kussmaul, 
2000, p. 83 y 84). 
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Figura 3: A) Contexto geológico regional del área de estudio dentro del 
Grupo Irazú. SFBH: Sistema fallas Bajo La Hondura, SFRS: Sistema 
fallas Río Sucio, FAG: Falla Alto Grande, FB: Falla Blanquito, FLi: Falla 
Liebres, Fl: Falla Irazú, FL: Falla Lara, FC: Falla Cipreses. B) Unidades 
litológicas locales identificadas dentro de la estructura volcánica Pico de 
Piedra. 
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Fuente: Alvarado et al. (2006), Bergoeing (1979), Hidalgo, Alvarado & Linkimer (2004), 
ITCR (2014), Linkimer (2003), Montero & Alvarado (1995), Montero et al. (2013), Ro- 
dríguez, Linkimer & Montero (2019) y Vargas (1994) 


Tectónicamente, el sitio es muy dinámico y presenta varias fallas 
activas a las cuales se les ha atribuido sismicidad histórica e instrumental 
(Montero & Alvarado, 1995, p. 28 y 29; Peraldo y Montero, 1999). Según 
Montero et al. (2013, p. 868) las fallas del flanco norte del volcán Irazú 
conforman el Cinturón Deformado del Centro de Costa Rica y han causa- 
do sismos destructivos atribuidos al Sistema de Fallas Río Sucio y Lara, 
respectivamente (Figura 3). 
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Estos sismos han alcanzado intensidades Mercalli Modificadas de 
hasta IX (Montero & Alvarado, 1995, p. 30; Peraldo y Montero, 1999, p. 
228) o aceleraciones espectrales estimadas de hasta 2,1 m/s? (21% de g) 
según modelos empíricos de intensidad (Linkimer, 2008, p. 84). 

Desde la amenaza símica, el registro más lamentable es el impacto 
que generó el Terremoto de Patillos del 30 de diciembre de 1952 en el Sis- 
tema de Fallas Río Sucio, causando 21 muertes, un impacto geológico y 
geotécnico amplio en el área mesosísmica, como deslizamientos, inestabi- 
lidad en estructuras y flujos de lodo, impactando localidades como Casca- 
jal, El Naranjo, Tierras Morenas, Rancho Redondo y San Isidro (Montero 
& Alvarado, 1995, p. 31; La Nación, 1953, p. 6). Otro evento es el Terre- 
moto de San Antolín, Cartago de 1841 que generó 22 muertes, destruc- 
ción masiva en Cartago y deslizamientos importantes en las localidades 
de Vázquez de Coronado y Goicoechea, asociándose según descripciones 
históricas con la Falla Lara (Peraldo y Montero, 1999, p. 227). 

Por tanto, son las fallas Río Sucio y Lara principalmente, las que 
pueden generar sismos que disparen inestabilidad de laderas, caída de ro- 
cas o flujos de lodo como se han reportado históricamente para el cauce del 
río Virilla (Alvarado, 2011, p. 257). 


Metodología 

Considerando el contexto geológico y de amenaza sísmica e hidroló- 
gica en la cuenca alta del río Virilla, se enfatizó en caracterizar los proce- 
sos de geodinámica externa, entre ellos la ocurrencia de deslizamientos e 
inventariado de deslizamientos actuales. 

Para esto se empleó teledetección, vistas al campo, caracterización 
geomorfológica del área, así como el modelado de susceptibilidad y esta- 
bilidad de laderas en la cuenca a partir de la metodología Mora-Vahrson 
modificada (Mora, Vahrson & Mora, 1992) y Probabilistic Infinite Slope 
Analysis (PISA) propuesta por Haneberg (2005). Los insumos geográficos 
empleados fueron la cartografía a escala 1:10.000 del Instituto Geográfi- 
co Nacional (IGN), las ortofotos del Sistema de Información Territorial 
(SNIT) a escala 1: 5.000 de 2017. El manejo e interpretación de insumos y 
datos espaciales se realizaron en el software ArcGIS y QGis. 

La escogencia de estas metodologías para estimación de ocurren- 
cia de deslizamientos, permite comparar las áreas susceptibles y con 
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probabilidad de afectación, ya que en el caso de la metodología Mora- 
Vahrson Modificada (MVm), la asignación de pesos y sus parámetros se 
asignan con categorizaciones cualitativas a criterio de experto y por su na- 
turaleza no considera elementos de probabilidad ni aproxima algún factor 
de seguridad al sector en estudio (p.e. Mora, Chávez & Vázquez, 2002). 

La metodología PISA permite incorporar parámetros geotécnicos 
de materiales y considera un coeficiente sísmico de diseño, por lo que de- 
termina en condiciones estáticas y pseudoestáticas el factor de seguridad 
de la ladera por unidad de área (pixel) y la probabilidad de excedencia de 
aceleraciones críticas del suelo; sin embargo, posee algunas limitaciones 
que se discutirán posteriormente. 

El pixel en este caso tiene una resolución de 1,5 m en dependencia de 
la información topográfica disponible, por lo que los parámetros restantes 
de las metodologías MVm y PISA se ajustan al menor tamaño de celda 
posible, inclusive las estimaciones de intensidad sísmica e hidrológica. 

La aplicación de ambas metodologías en la cuenca alta de río Viri- 
lla, constituye una oportunidad comparativa muy particular para estimar 
la afectación por procesos de remoción en masa en un ambiente cerrado y 
controlado, como lo es la estructura volcánica Pico de Piedra. 

La metodología Mora-Vahrson modificada (Mora, Vahrson & Mora, 
1992) es un método heurístico que considera a priori cinco elementos pa- 
rámetros para la estimación de susceptibilidad al deslizamiento, valién- 
dose de mapas geológicos, humedad del suelo, pendiente del terreno y 
disparo por sismos o lluvias intensas para un periodo de retorno de 100 
años. Por tanto, requiere de instrumentación meteorológica cercana, reg- 
istros sísmicos históricos y caracteriza el grado de susceptibilidad (H) en 
muy alta, alta, moderada, baja y muy baja (Fórmula 1). 


H = (S; * Sı* Sn) * (T; +T,) [1] 


Donde: 

Ss: Susceptibilidad por pendiente 

Sl: Susceptibilidad litológica 

Sh: Susceptibilidad por humedad del suelo anual 

Ts: Disparo por sismo 

Tr: Disparo por lluvia (Duración: 60 min Tretorno= 100 años) 
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Las pendientes se obtuvieron a partir de datos topográficos del IGN, 
mientras que datos de humedad del suelo y disparo por lluvia se obtuvie- 
ron de registros abiertos del CATIE y del Instituto Meteorológico Nacional 
(IMN). Algunos casos exitosos de aplicación en Costa Rica se presentan 
en Segura, Badilla & Obando (2011), Mora, Chávez & Vázquez (2002), 
Barrantes, Barrantes £ Núñez (2011) y Quesada & Feoli (018, p. 34). 
Para el detalle de puntuación por parámetro puede consultarse Segura, Ba- 
dilla ££ Obando (2010, p. 106-114); respecto a los parámetros empleados 
en esta investigación se especifican en el Cuadro 1, considerando inten- 
sidades sísmicas similares a las del Terremoto de Patillos de 1952. 


Cuadro 1. Parámetros empleados del método MVm en las unidades 
litológicas identificadas. 


Parámetros de entrada Descripción van 
asignado 
Susceptibilidad por pen- o 5 
diente (Ss) De 0°a71 0a6 


Suelos volcánicos regolíticos, rocas volcánicas 


Susceptibilidad litológica andesíticas y piroclásticas, coluviones y depósitos 3,4y5 


(8) aluviales. Nivel freático alto (Figura 3 B). 
Susceptibilidad humedad Anual, estaciones cercanas IMN. 3 
suelo (S) De 2325 a 2630 mm. 
Disparo por sismo (Ts) IMM = VIII, aceleración pico = 0,21% g. 6 
Duración 60 min, T somo = 100 años. 
Disparo por lluvia (Tr) Estación CATIE (1949-2016). 4 
Distribución Gumbel,P_.. = 368 mm/hora. 


máxima 


Fuente: Mora, Chávez & Vázquez (2002), Sequeira-Arguedas (2019), IMN (2019) y Montero & 
Alvarado (1995) 

La metodología Probabilistic Infinite Slope Analysis (PISA) considera 
para sus cálculos la ecuación del talud infinito (Haneberg, 2005, p. 14), por 
tanto solo considera espesores deslizables de suelo y rocas como láminas o 
fallas planares (González de Vallejo, Ferrer, Ortuño & Oteo, 2002, p. 440 y 
447-449). Es una metodología predictiva que presenta fundamento geomor- 
fológico, geotécnico y sobretodo estadístico para la asignación del factor de 
seguridad en condiciones estáticas o la probabilidad de excedencia de las 
aceleraciones críticas de Newmark (el suelo alcanza un factor de seguridad 
< 1,0) en caso de considerarse un sismo de diseño (Fórmulas 2, 3 y 4). 
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De acuerdo con la distribución estadística de cada variable aleatoria 
o parámetro que conforman la ecuación del talud infinito (Fórmula 2), se 
emplean aproximaciones e iteraciones generadas con estimación inferen- 
cial puntual (First Order Second Moment) de los estadísticos de la muestra 
(parámetros), aproximando una iteración de Montecarlo para un estadísti- 
co buscado (media aritmética del FS). 

Se genera un archivo espacial ráster (grid file) de estabilidad en tér- 
minos del factor de seguridad (FS) o probabilidad de excedencia de una 
aceleración pico dada (PGA) bajo condiciones estáticas y pseudoestáticas 
(Fórmulas 3 y 4) para cada unidad geotécnica del área, el cual es desplega- 
ble en cualquier sistema de información geográfica (SIG). Este análisis no 
considera falla de la ladera por deslizamientos rotacionales o mecanismos 
compuestos y la Figura 4 ejemplifica el elemento de suelo o roca desliza- 
ble que considera la metodología PISA. 


Eg = ECH iD Hat Yin By PICOS P tanp. 
[qt + YmD+(Ysat T Ym) Hw D] cos B sin fP 


crit = (FS — 1)g sin B [3] 


[2] 


Play < acrit] = F [ay (acrit)] [4] 


Las ventajas de este método, es que permite incorporar parámetros 
geotécnicos, asignar carga y cohesión de la vegetación, además de que 
puede modelarse cualquier escenario de saturación del suelo (seco, semi- 
saturado y saturado) y cualquier evento sismico (% de g) a partir de un 
coeficiente de aceleración pico o sostenida. Algunos casos exitosos de 
aplicación de la metodología PISA son Mora (2010, p. 48 y 49), Granados 
(2013), Quirós (2019) y Carrillo (2019, p. 104). A nivel internacional se 
tienen investigaciones como Griffits, Huang & Fenton (2011). 

El rango de variación por parámetro se recopiló en literatura a nivel 
nacional e internacional, empleando una media y su desviación estándar 
para cada unidad geotécnica (Cuadro 2). Las unidades empleadas son las 
establecidas por el Sistema Internacional de Unidades (SI). 
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Figura 4. Elemento de suelo o roca deslizable considerado en la ecuación 
del talud infinito. 


Parámetros empleados 


C: cohesión aportada por raíces (Pascal, Pa) 
C? cohesión del suelo (Pa) 

qt carga aportada por la vegetación (Pa) 
Ym peso unitario húmedo del suelo N/m) 
Ysat peso unitario saturado del suelo Nin ) 
Yw peso unitario del agua (9810 Nim ) 

D espesor de suelo deslizable (m) 
Hw coeficiente saturación suelo (De 0 a 1) 
15 ángulo de la ladera (?) 

p' ángulo de fricción interna del suelo (°) 


AN ACrit aceleración inicial y crítica de Newmark (m/s) 
Fuente: González de Vallejo et al. (2002) y Haneberg (2005) 


Cuadro 2. Parámetros empleados del método PISA para las unidades 
geotécnicas identificadas. 


Distribución 


Parámetro Rangos considerados a Referencia 
estadística 
Entre 0 y 3000 Pa (bosque pe- 
C, (Pa) renne y coníferas en suelos de Uniforme López et al. (2008) 
montaña) 
De 0 a 0,066 Pa (57% área de 
q, (Pa) bosques) Media: 0,0335 Pa Triangular López et al. (2008) 
Arboles de 4 a 6 m altura 
13 700 + 2400 (Unidad SV) Mora (1998), Granados 
Y, (KN /m>) 19 120 + 3220 (Coluviones) Noral (2013), Jiménez et al. 
18 850 + 3220 (Aluvión) (2014, 2014b) 
15 600 + 1700 (Unidad SV) Mora (1998), Granados 
Ya (KN/m’°) 20 340 + 2900 (Coluviones) Normal (2013), Jiménez et al. 
20 140 + 2900 (Aluvión) (2014°, 2014b)* 
y, N/m’) 9810 N / m? Constante - 
D Media: 7+5m Normal Inspección en campo 
H, De 0 a 100 % (0a 1) Uniforme e 
p Curvas de nivel cada 5 m Variable Curvas de nivel I.G.N. 
6000 + 6000 (Unidad SV) Mora (1998), Granados 
C (Pa) 113 700 + 113 700 (Coluviones) Normal (2013), Jiménez et al. 
29 400° + 20 000 (Aluvión) (Q014*, 2014b)* 
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Distribución 


Parámetro Rangos considerados E Referencia 
estadística 
33 + 6,7” (Unidad SV) Mora (1998), Granados 
p’ 25.7° + 10° (Coluviones) Normal (2013), Jiménez et al. 
26,7° + 9,71° (Aluvión) (2014, 2014b)* 
> Datos LIS, Montero & 
a, Y a o 8 Constante Alvarado (1995), Lin- 


kimer (2008) 


Fuente: Unidades geotécnicas basadas en Sequeira-Arguedas (2019) y Mora (1998) 


Respecto a la evaluación ante flujos de lodo y avenidas extraordina- 
rias, se recorrió el cauce del río Virilla para las secciones cercanas al corredor 
vial La Holanda y se evaluó la granulometría de los depósitos, discriminan- 
do una zona de influencia a partir de la morfología del cauce y alcance de 
depósitos de inundación. No se realizó un análisis hidrológico exhaustivo. 


Resultados y discusión 

Geomorfología 

La cuenca en estudio presenta un perímetro de 7,97 km, 2,88 km? de 
área, largo de 2,96 km y un ancho máximo de 1,62 km. El sitio presenta 
morfologías volcánicas con predominancia de pendientes altas (en algunos 
sectores hasta 70%), además de evidencias de inestabilidad como repta- 
ción, deslizamientos superficiales, deslizamientos antiguos y afloramiento 
de depósitos coluviales en inspecciones de campo, las cuales se clasifican 
como formas de origen denudacional (Figura 5 y 6). 

Además se identificaron, formas fluviales asociadas a la dinámica 
del río Virilla y sus tributarios, generando pequeñas áreas de deposición 
fluvial y arrastre de bloques rocosos. 

La evolución de la estructura volcánica Pico de Piedra generó las 
morfologías de escarpes, colinas volcánicas de moderada y alta pendiente, 
además del cono intracratérico secundario, localizado en el extremo sures- 
te del área de estudio (Fig. 5). 

Estas unidades presentan patrones de drenaje paralelos y diversos 
grados de disección, alcanzando profundidades de drenajes de hasta 70 m, 
aunque la disección es mayor en los valles externos a la cuenca del Virilla. 

Las colinas con pendientes superiores a 35° usualmente presentan 
reptación, divisorias agudas y deslizamientos superficiales recientes, ade- 
más en época lluviosa se activan cauces ocasionales que arrastran material 
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hacia los valles fluviales o son depositados en las zonas de piedemonte o 
de transición (Figura 6). La erosión es intensa y deja algunos afloramientos 
del macizo rocoso, variando el perfil de suelo entre 1 y 12 m de espesor. 

Las terrazas intracratéricas son remanentes de coladas de lava y piro- 
clastos de la Formación Reventado, que por su grado de alteración y mor- 
fografía pueden asociarse a uno o varios eventos volcánicos, presentando 
pendientes promedio de 15° en las secciones de rellanos y pendientes de 45° 
o más en los escarpes. Presentan reptación y son empleadas para ganadería. 

Las formas denudacionales son las modeladoras del relieve actual en 
la cuenca alta del río Virilla, ya que se idenficaron al menos tres desliza- 
mientos antiguos con coronas de 200 m de ancho y depósitos remanentes, 
además de nueve deslizamientos activos de anchos inferiores a 60 m con- 
centrados en las cercanías de los valles fluviales. 

Adicionalmente, se identifica un cerro disectado en el sector central 
de la cuenca cuyo modelado se atribuye a agentes hídricos que aceleraron 
la meteorización química en las rocas, aunque su génesis también puede 
asociarse a colapsos gravitacionales que aceleraron el proceso de disección. 

Los depósitos de piedemonte se localizan en las zonas deprimidas 
cercanas a los escarpes volcánicos, presentando ondulaciones tipo hum- 
mocky, bloques angulosos métricos inmersos en suelo, alcanzando el co- 
rredor vial La Holanda y el drenaje de la cuenca (Figura 6). 

En el caso de las morfologías fluviales, se restringen a los talweg del 
río Virilla. En los sectores de mayor elevación se presentan valles en “V” 
con secciones de hasta 25 m de ancho, mientras que en los sectores más 
bajos predominan valles fluviales con secciones transversales de hasta 60 
m de ancho y sectores inundables. En este sector se reconocieron entre uno 
y tres niveles de terrazas fluviales apareadas las cuales presentan alturas 
entre 75 cm y 1,5 m (Figura 7). Los depósitos presentan imbricación y 
soporte por matriz arenosa. 

Sin embargo, se localizaron bloques erráticos y vegetación arrastra- 
da en los rellanos de las terrazas (Figura 7), por lo que estos sectores son 
muy dinámicos y dependen de la escorrentía en la cuenca, particularmente 
en los sectores más elevados. 
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Figura 5. Mapa de unidades geomorfológicas identificadas para el 
área de estudio. Los triángulos amarillos denotan facetas triangulares 
estructurales inferidas. 


506000 507000 508000 


- Formas de origen volcánico 
l-- -i a [C] Colinas a disectadas alta pendiente (>35°) EL Deslizamiento antiguo 

=— Falla geológica activa Colinas volcánicas moderada pendiente (10- 359) al Depósito deslizamiento antiguo 
Ríos y quebradas E Cono intracratérico residual Do e Piedemonte 

~ Sitios con reptación Formas de origen fluvial 

MM Terrazas fluviales río Virilla 
E] valles en “v” 


TJ] valle fluvial río Virilla y tributarios 


==] Escarpe interno de caldera residual 
A) Deslizamientos 


-] Terrazas intracratéricas de colapso (x3) 


* SRC: CRTM-05 p P z 
ua CRUS Formas de origen denudacional 


Elipsoide: WGS-84 [EA] Cerro residual volcánico 


Fuente: Elaboración propia 


Mora-Vahrson modificado 

La susceptibilidad ante deslizamientos en la cuenca alta del río Vi- 
rilla según la metodología Mora-Vahrson modificada permite evaluar de 
manera relativa los sectores más inestables. La inestabilidad es más ex- 
tendida en los sectores de mayor pendiente si se considera un sismo con 
intensidades Mercalli Modificadas de VIII (Figura 8), ya que bajo una pre- 
cipitación atípica de hasta 368 mm en una hora, la mayoría del área se cla- 
sifica como de susceptibilidad moderada y escasos sectores como de baja, 
relacionado al efecto de la pendiente. 
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Figura 6. A) y B) Vista del escarpe volcánico, las zonas de deposición 
coluvial y cauces ocasionales. C) Escarpe volcánico de 150 m y coronas 
de deslizamiento observadas. 


Zona de deposición 
coluvial 


- Zona de deposición l 
j > E, coluyial > es 0% X Càuce fio Virilla 


Fuente: Elaboración propia 


Figura 7. A) Niveles de terrazas aluviales en el margen izquierdo del 
río Virilla (aguas abajo). B) Vista de los materiales arrastrados por el 
río. C) Vista de la granulometría típica de depósitos aluviales en niveles 
antiguos. D) Vista de niveles más antiguos de terrazas. 


na alcance depósitoliyales 


Bloques rocosos y 
troncos en el cauce 


Fuente: Elaboración propia 
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Existe un escenario crítico si sucede un sismo cercano como el even- 
to de Patillos de 1952 y si los suelos se encuentran saturados (p.e. entre 
septiembre a diciembre), ya que en los escarpes y terrazas volcánicas la 
susceptibilidad es muy alta. Otros sectores en la cuenca se clasifican como 
de alta y moderada susceptibilidad en función de los parámetros del méto- 
do (Fig. 8), siendo susceptibles las áreas de deslizamientos antiguos. 

Un sismo localizado en la Falla Lara sería un escenario crítico, pues 
su cercanía podría generar aceleraciones en el suelo superiores a las fuer- 
zas estabilizadoras intrínsecas de los materiales, ya que este escenario es 
posible considerando la Falla Lara como el foco del sismo de San Antolín, 
Cartago de 1841 (Peraldo & Montero, 1999, p. 227 y 228), el cual generó 
deslizamientos considerables en la región. 


Figura 8. Susceptibilidad al deslizamiento en la cuenca alta del río 
Virilla, según la metodología Mora-Vahrson modificada. 


Susceptibilidad (H) 
E] Muy baja (<51) Deslizamientos 
E] Baja (51-142) —— Calles 
[L _] Moderada (143-288) 
[UN] Alta (289-430) 
EE] Muy alta (> 430) 


— Ríos 


Fuente: Sequeira-Arguedas (2019) 


Probabilistic Infinite Slope Analysis (PISA) 

Considerando los modelos de estabilidad y con las propiedades 
geotécnicas y físico-mecánicas de las unidades geotécnicas (Fig. 3B): A) 
suelos volcánicos poco a moderadamente compactos (SV), B) aluviones 
con matriz arenosa, C) Coluviones y material piroclástico suelto; (Cuadro 
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2); se decidió emplear un escenario en condiciones sin sismo (estáticas), 
un escenario con suelo saturado, otro con un sismo lejano en el arco vol- 
cánico similar al evento de Cinchona de 2009 y otro con un sismo cercano 
con aceleraciones críticas de hasta 2,06 m/s? (206 Gal) con similitudes al 
sismo de Patillos. 


Figura 9. Media del factor de seguridad en las laderas en condiciones de 
suelo seco y saturado y su grado de estabilidad asociado, empleando la 
metodología PISA. 


Condición estática 
Media del Factor de 
Seguridad Estático ladera (FS) 


Hu < 1,0 (Muy Inestable) 
ES 1.0- 1.2 (Inestable) 
EST 1.2 -1,4 (Mod. Estable) 
MA  —>1.4 (Estable) 


f- + 
3 ø ha 
Suelo seco 
Sin sismo 


% del área 


E 
pu 


"i fa P. e 
Suelo saturado X 
f Sin sismo ~ 


E E i - 


Seco Saturado 


Fuente: Sequeira-Arguedas (2019) 


En condiciones estáticas (Figura 9) es notorio que las zonas inesta- 
bles se restringen a los escarpes (23% de área), que por su condición geo- 
morfológica, poseen pendientes mayores al ángulo de fricción interno de 
los materiales, manteniendo una similitud con el caso de disparo por sismo 
de la metodología Mora-Vahrson modificada. 

El caso crítico ocurre en condiciones de suelos saturados, la cual es 
probable durante los meses de invierno en la cuenca (septiembre-diciem- 
bre) y si existe una concentración anómala de precipitación por efecto de 
una tormenta. En este escenario, solamente un 12% del área puede catalo- 
garse como estable, ya que más de un 75% de la cuenca alcanza factores 
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de seguridad cercanos a 1,0, pudiendo generar deslizamientos y aumento 
en los procesos erosivos ya identificados. 


Figura 10. Probabilidad de excedencia de la aceleración crítica en las 
laderas (FS < 1,0), en dos sismos o escenarios críticos empleando la 
metodología PISA. 


Condición pseudoestática 
Probabilidad excedencia (%) de 
la aceleración crítica del suelo 
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Fuente: Sequeira-Arguedas (2019) 


El caso pseudoestático (Figura 10) implica considerar una fuerza ex- 
terna adicional en las laderas, la cual generalmente se asume paralela al 
vector de máxima pendiente y cuya magnitud depende de las aceleraciones 
sísmicas espectrales en el sitio. A pesar de que cualquier escenario sismico 
puede modelarse, se consideran dos eventos históricos generados en el 
arco volcánico de la Cordillera Central. 

En el caso del sismo de Cinchona de 2009, considerando una acele- 
ración de 0,98 m/s? solamente las áreas de alta pendiente son las que alcan- 
zarían eventualmente inestabilidad, debido al efecto de atenuación sísmica. 
Cerca de 21% del área presenta una probabilidad alta de alcanzar un factor de 
seguridad pseudoestático inferior o igual a 1,0. Se modelaron las condicio- 
nes de saturación aproximadas en enero según los registros de precipitación; 
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sin embargo, al tratarse de suelos volcánicos estos podrían encontrarse con 
contenidos de humedad mayores o casi saturados (Hw ~ 1,0). 

En el caso de un sismo cercano con suelo saturado como el de Pa- 
tillos en diciembre de 1952, el escenario de estabilidad es opuesto, pues 
cerca de un 82% del área de la cuenca tiene una probabilidad alta a muy 
alta de alcanzar la inestabilidad. Esta probabilidad indica que un pixel (de 
1,5 m) puede tener una probabilidad del 50% o mayor de deslizarse, según 
el modelo de falla escogido (talud infinito). Solamente las áreas corres- 
pondientes a las terrazas volcánicas y algunos valles fluviales pueden con- 
siderarse “estables”; sin embargo, en un análisis integrado por amenaza 
geológica no es viable considerar estos sitios como poco susceptibles, pues 
los cauces pueden represarse producto de un sismo intenso o ser propensos 
a caída de materiales, entre ellos bloques, árboles o suelo inconsolidado. 
Registros históricos evidencian que la cuenca del río Blanco (4 km al norte 
del área de estudio) se vio afectada por represamientos y avalanchas du- 
rante el sismo de Patillos de 1952 (Montero & Alvarado, 1995, p. 29-35), 
por lo que es posible un escenario similar en el área. 


Capacidad de arrastre del cauce y su relación con los puentes del 

corredor vial 

En visitas al campo, se comprobaron bloques rocosos ígneos angulo- 
sos y subredondeados, que generan alta rugosidad del cauce, presencia de 
bancos laterales de deposición y sectores con incisión del cauce en mate- 
riales blandos, principalmente arcillas (Figuras 11 y 13). 

Inspeccionando los tramos superiores e inferiores a los puentes 1 y 
2, se encontraron bloques rocosos cuyas dimensiones superan los claros 
actuales o altura libre respecto al espejo de agua actual, con lo cual se 
determina que flujos concentrados de detritos han ocurrido con anteriori- 
dad y pueden generarse nuevamente, siendo capaces de movilizar enormes 
volúmenes de material y dañar seriamente los accesos al sitio, siendo una 
zona estratégica de producción de lácteos y hortalizas, además de captarse 
agua para consumo humano. 

Durante la actividad eruptiva del volcán Irazú entre 1963 y 1965, por 
el cauce del río Virilla descendieron lahares que alcanzaron el poblado de 
San Isidro de Coronado y otros poblados aguas abajo en el mes de septiem- 
bre de 1963 (Alvarado, 2011, p. 257). Las subestructuras de los puentes 
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son susceptibles por impactos, arrastre hidráulico y oxidación en sus com- 
ponentes, principalmente el Puente 2, que presenta un estado estructural 
pobre, con varillas expuestas, apoyos agrietados y vigas corroídas. 

El Puente 1 presenta muy buen estado de servicio, pues su estructura 
de concreto reforzado no presenta agrietamientos ni alteración. La antigúe- 
dad de estos puentes es de 50 años y ya han soportado eventos sísmicos 
(Cuadro 3). 


Figura 11. A) Puente 1 en el corredor vial La Holanda, fabricado en 
concreto reforzado. B) Puente 2 en el corredor vial La Holanda, fabricado 
en acero, madera y subestructura de concreto reforzado. C) Bloque 
rocoso anguloso en márgenes del cauce (escala 1 m). 


MA Tendido 
eléctrico 


ma 


y 
Subestructuta. 7 


Fuente: Elaboración propia 
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Cuadro 3. Especificaciones de los puentes 1 y 2 en el corredor vial La 


Holanda. 
Lar Anch 1 ; 
Puente Estructura g0 Se S PA Observaciones 
(m) (m) Libre (m) 
Superestructura Cimentado en depósi- 
ti a i = 
Puente 1 y subestructu 15,7 58 28 tos aluviales, obstruc 
ra de concreto ciones de vegetación. 
reforzado. Buen estado. 
: F a 
Superestructura de Cimen ado ewdepó 
sitos aluviales. Acero 
agoro y madera. corroído. Presenta 
Puente 2 Subestructura 11,8 3,1 2,5 ON ] 
filtraciones y agrieta- 
de concreto . 
mientos. Pobre estado 
reforzado. s5 
de servicio. 


Fuente: Elaboración propia 


Zonificación y recomendaciones de uso 

Considerando la capacidad de arrastre hidráulica, la amenaza sís- 
mica, la susceptibilidad y ocurrencia probabilística de deslizamientos en 
la cuenca alta del río Virilla, se generó una zonificación integrada (Figura 
12), con el fin de determinar cuáles sectores presentarían afectaciones si 
ocurre un evento sísmico o hidrometeorológico. 

En primera instancia se generó una zona inundable (retiro) de todos 
los cauces de la cuenca con 50 m de ancho y se determinaron las áreas de 
alcance por colapso gravitacional y deposición coluvial. 

Con la susceptibilidad y estimación del factor de seguridad en las 
laderas, se discriminaron las áreas sumamente inestables, moderadamen- 
te estables y de uso conforme. Las áreas inestables coinciden con las 
de mayor pendiente y escarpes con reptación, en ellas es recomendable 
conservar las zonas forestales con el fin de disminuir el potencial erosivo 
por escorrentía. 
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Figura 12. Mapa de zonificación integrado ante amenazas geológicas 
y uso recomendado en la cuenca alta del río Virilla, considerando 
escenarios sísmicos críticos y suelo saturado. 
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Las áreas moderadamente estables son escasas y actualmente se de- 
dican a actividades ganaderas; sin embargo, no se recomienda remover 
zonas forestales, sino más bien controlar la erosión fluvial, de ser posible 
instrumentar la ladera en lo que respecta a las haciendas productivas. En 
las áreas de uso conforme se ubica infraestructura habitacional y producti- 
va, representando un 24,6% del área. 

Estos sectores son estables pero pueden presentar algún grado de 
afectación ante caída de materiales. En el caso del corredor vial La Ho- 
landa, se determinó que 2,1 km son susceptibles a colapsar (un 39,6% del 
corredor vial), pues presentan socavaciones laterales producto de la acción 
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erosiva del río Virilla, depósitos de caída, desprendimientos activos y es- 
casos drenajes pluviales (Figura 13). 


Figura 13. A), B) y C) Vista de desprendimientos recientes a lo largo 
del corredor vial La Holanda que atraviesa el área de estudio, notorio los 
escasos desagúes pluviales. D) Claro vertical de 2,8 m en el Puente 1 y 
materiales represados. 
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Fuente: Elaboración propia 


Conclusiones 

La cuenca alta del río Virilla es el resultado de la evolución de un 
foco eruptivo antiguo, asociado al volcanismo del Grupo Irazú, cuyas rocas 
datan del Chibaniano o Pleistoceno Medio (Alvarado et al., 2006, p. 263 y 
264) y el modelado por la acción del río Virilla y eventos de colapso acti- 
vos e inactivos. Predominan laderas de pendientes superiores a 10°, en sec- 
tores escarpados donde la erosión es intensa y la inestabilidad es evidente, 
con presencia de reptación y deposición coluvial en los piedemontes. Esta 
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serie de agentes de geodinámica externa ha generado disección, divisorias 
sinuosas agudas y patrones de drenaje paralelos. 

Las metodologías Mora-Vahrson modificada (MVm) y PISA cons- 
tituyeron un insumo importante en la determinación de la susceptibilidad 
al deslizamiento en la cuenca así como evaluación probabilística de esta- 
bilidad en la ladera, siendo posible complementarlas y modelar diversos 
escenarios sísmicos y de saturación de suelos. Sin embargo, MVm es útil 
en estudios a priori y es necesario contar con un modelo geológico detalla- 
do. PISA por su parte, trabaja en términos probabilísticos pero requiere de 
una asignación detallada de unidades geotécnicas (espesores deslizables y 
parámetros de corte), con la ventaja de que permite incorporar propiedades 
de resistencia al corte y físicas de los materiales. 

Con la integración de factores geológicos evaluados, existe con- 
gruencia con el uso actual y las zonas de mayor estabilidad en la cuenca; 
sin embargo, no implica que no se verán afectadas a futuro ante sismos 
cercanos y eventos hidrometeorológicos extraordinarios. Un tramo impor- 
tante (2,1 km) del corredor vial La Holanda es susceptible a colapsar por 
desprendimientos de material, flujos de lodo (puentes) y escasos desfogues 
pluviales. Es necesario conservar la cobertura forestal para disminuir el 
potencial erosivo por escorrentía. 
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